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Abstract 

We discuss whether or not "freezing" of the QCD running cou- 
pling constant in the infrared region is consistent with the Schwinger - 
Dyson (SD) equations. Since the consistency of the "freezing" was not 
found, the conclusion is made that the " analytization" method does 
not catch an essential part of nonperturbative contributions. Proceed- 
ing from the results on consistency of the infrared enhanced behaviour 
of the gluon propagator with SD equations, the running coupling con- 
stant is modified taking into account the minimality principle for the 
nonperturbative contributions in the ultraviolet region and conver- 
gence condition for the gluon condensate. It is shown that the re- 
quirements of asymptotic freedom, analyticity, confinement and the 
value of the gluon condensate are compatible in the framework of 
our approach. Possibilities to find an agreement of the enhanced be- 
haviour of the running coupling constant with integral estimations in 
the infrared region are also discussed. 



^Talk presented at the Xllth International Workshop on High Energy Physics and 
Quantum Field Theory, September 4 - 10, 1997, Samara, Russia 



The phenomenon of asymptotic freedom |I| called forth an impressive suc- 
cess of perturbative QCD in the description of experimental data plethora. 
However, there is wide scope of phenomena which is intractable in the frame- 
work of perturbation theory. Nonperturbative effects modify the infrared 
behaviour of the quark and gluon Green's functions. With decreas the 
renormalization group improved one-loop running coupling constant 

(60 = IIC2/3 — 2A^//3) increases, which may indicate a tendency of unlimited 
growth of the interaction at large distances, leading to a confinement of 
coloured objects. However, at = Is? in the pole is present, which 
is nonphysical, at least, due to the fail of the perturbation theory, and the 
account of nonperturbative effects becomes obligatory. 

In recent papers the solution of the problem of a ghost pole was pro- 
posed with the condition of analyticity in q^ being imposed. The idea of 
"forced analyticity" goes back to p, H of the late 50s. They were dedicated 
to the problem of Landau-Pomeranchuk pole in QED. Using for 6is{q^) 
a spectral representation without subtractions, the following expression for 
the running coupling constant was obtained in O] 



+ 



In(g2/A2) A2 - q 



(2) 



This expression is analytic in the infrared region due to nonperturbative 
contributions and it has a finite limit at zero (although the derivative is 
infinite) . 

Nowadays a possibility of "freezing" the coupling constant at low energies 
is under discussion in the framework of some scheme of approximate 
calculations ((16| — Nf) - expansion). In the approach [0 with confining 
background field, "freezing" was also obtained, 

47r 

"^^^'^ ^ 6oln((m| + g2)/A2)- (3) 

Here is a process-dependent constant of the order of 1 GeV. 

In the present paper we follow the approach of Refs. P, P]. We discuss 
the problem of consistency of the constant behaviour of the running coupling 
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constant in the infrared region with Schwinger- Dyson (SD) equation for a 
gluon propagator. Further we include into consideration additional nonper- 
turbative terms, in particular, the singular in the infrared region term ~ 
the necessity of renormalization invariance being taken into account. Then 
we consider possibilities to fulfill the demands of confinement, asymptotic 
freedom, analyticity, correspondence with estimates of the gluon condensate 
and integral estimates for Q;s(g^) in the infrared region. 

Let us consider the integral SD equation for the gluon propagator in 
ghost-free axial gauge [|10] A^r]^ = 0, 77^ — gauge vector, rf 7^ 0. In this 
gauge the running coupling constant is directly connected with the gluon 
propagator and Slavnov- Taylor identities have the simplest form. An 
important preference of the axial gauge consists in a possibility to exclude 
the term from the SD equation, which contains a full four-gluon vertex, 
by means of contraction of the equation with tensor ri^rii,/rf. We work in 
the Euclidean momentum space, where smallness of the momentum squared 
means smallness of all its components. For the gluon propagator Dfj_^{p), we 
suppose the approximation D^y{p) = Z{p'^)D(^Q)^^{p) to be appropriate for 
studying the infrared region. D(^Q)^^{p) is a free gluon propagator. In the first 
paper of Ref. the possibility of the infrared behaviour for renormalized 
function = Z{p'^) / Z{^'^) being of the form 

Zr{x)=Zr{{])+o{1), a:^+0, (4) 

has been studied {Zf{{0) is nonzero constant) in the framework of nonpertur- 
bative approach of Baker - Ball - Zachariasen (BBZ) [0. The approxima- 
tions of the BBZ approach as well as the condition y = imposed on gauge 
parameter y = {priY /p^rf seems to be adequate to studying the possibility 
of the infrared behaviour which is not too singular. With the assumption (^) 
the equation for function Z^(p^) takes the form 



Z-R^p") = 1 + Zniof-^^ W + o{\np'),p' 0. (5) 



.^11 
167r2 3 

We see, that the behaviour ZRijP') ^ Z{Q) 7^ for — > does not agree 
with the SD equation. 

This conclusion encourages us to look for the possibilities different from 
the assumption on the finiteness of the coupling constant at zero. Recently 
a possibility of the soft singular power infrared behaviour of the gluon prop- 



agator has been discussed |T3[, D{q) ~ (g ) q — >• 0, where /3 is a small 
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positive non-integer number. In Ref. [|T^] the consistency of such behaviour 
with the same equation was studied. A characteristic equation for the expo- 
nent (3 was obtained and this equation was shown to have no solutions in the 
region < /5 < 1 . The authors of Ref. |jl5[ also came to the conclusion on the 
inconsistency of the soft singular infrared behaviour of the gluon propagator. 
The case of possible interference of power terms was studied in Ref. [1^ and 
it was shown that in a rather wide interval — 1 < /3 < 3 of the non- integer 
values of the exponent the characteristic equation had no solutions 0. At 
present a more singular, in comparison with free case, infrared behaviour of 
the form D{q) ~ MV(g^)^ ^ seems to be most justified [|17|, |T2|, |1^]. 
The models based on this assumption are widespread and rather attractive. 
In this way it is possible to describe the dynamical chiral symmetry breaking, 
to solve the U{1) problem, to evaluate the topological susceptibility |19|, pO|] , 



to calculate the condensates of gluon and quark fields |2T[], etc. The phys- 
ical consequences of such enhancement of zero modes are discussed in re- 
views p2| , ^ . Bearing in mind the remarks stated above, let us turn to the 



problem of nonperturbative contributions. The "analytized" expression (|^) 
without nonphysical singularity will be a starting point of our further con- 
sideration. We see that this expression has a nonperturbative tail with the 
behaviour at 0. To answer the question if this behaviour is ad- 

missible, let us consider the important physical quantity, namely, the gluon 
condensate K = < vac | Os/vr : G^j, G^j, :| vac > . According to the defini- 
tion (see e.g., [^) up to the quadratic approximation in the gluon fields, one 
has after the Wick rotation 

A- = - — y (^(■5..*' - KK)Dl^{k)yZ{k) - Z""(k)) = 

where a"^ is a nonperturbative part of the running coupling constant. The 
one-loop "analytized" behaviour of Eq. (||) leads to the quadratic diver- 
gence in Eq. (j^) at infinity and this is true for two- and three-loop expres- 
sions of the analytization approach. Note that at large g^, "freezed" 
behaviour (|^) and known enhanced in the infrared region Richardson's run- 
ning coupling constant (coinciding formally with (|^) at = A) do 

^Note that this interval contains values — 1 < /3 < for which the propagator vanishes 
at zero. 
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not ensure the convergence of the integral at infinity. According to the 
results of Refs. Jl^ |l^, |l^] let us add in Eq. (H) the isolated infrared singular 
term of the form This term is harmless at zero and it can improve the 

behaviour of the integrand at infinity and make the integral logarithmically 
divergent. To make integral convergent at infinity, it is sufficient to add 
one more isolated singular term of a pole type with parameters chosen ap- 
propriately. In this sense the model we come to is minimal. The expression, 
we obtain for the running coupling constant, is the following: 

47r / 1 A' cA2 (1 -c)An 

"'^^ ^ ~ 60 Vn(gVA2) ^ A2 -g2 + g2 + ^2 + ^2 j ' U 



with mass parameter m 
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m^ = AV(c-l). (8) 

It is worth noting that an account of nonperturbative contributions in Eq. (|^) 
preserves a perturbative time-like cut of Eq. (^. With the given value of the 
QCD scale parameter A, the parameter c can be fixed by the string tension 
K or the Regge slope a' = l/(27rK) assuming the linear confinement V{r) ~ 
Kr = a^r at r — > 00. We define the potential V{r) of static qq interaction 
[ p5| , p6| by means of three-dimensional Fourier transform of ots{(f')/(f' with 
the contributions of only one dressed gluon exchange taken into account. 
This gives the following relation 

cA^ = (36o/87r)a2 = {h^/lQ^^)g''M^ . (9) 

Taking a ~ 0A2GeV, one obtains c = A^/A^ where A^ = SboK/Sir ~ 
0.434: GeV {bo = 9 in the case of 3 light flavours). From Eq. (||) one ob- 
tains 

A2 

(10) 
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and the tachion absence condition limits the parameter A, A < ASAMeV. 
With Eq. d^) taken into account, the parameters in Eq. (|^ are only A and 
60- For SUcls) colour group {C2 = 3), iV/ = 3 and A = 300 MeV, the be- 
haviour of running coupling constant is represented in Fig. |I] for the cases 
under discussion. It should be noted that at sufficiently large q the curves 
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in Fig. |1] seem to be indistinguishable, but the presence of amphfying mul- 
tipler y in Eq. (H) allows one to distinguish the ultraviolet behaviour of the 
nonperturbative contributions. 

Let us represent expression (J^) in explicitly renormalization invariant 
form. It can be done without solving the differential renormalization group 
equations. In this order we write as{(f) = (p'isi^ll^-, Q^^I^t^ and use the 
normalization condition g^{\^g^^ = g'^. Then we obtain the equation for the 
wanted dependence of the parameter on g"^ and /i^: 

1 A2 A2 (l-c)A2' 

1 h C h 

ln(/i^/A^) A^ — /i^ /i^ + 

From dimensional reasons A^ = fi'^exp{—(f{x)}, where x = bog'^/lGir'^ = 
6o«s/47r, and for function ip{x) we obtain the equation: 

X = 1 h ce ' 

!f{x) 1 - e-^W (c - l)e'^(^) + 1 ■ 

The solution of this equation at c > 1 is function f{x), which has the be- 
haviour ip{x) ~ 1/x at X — i> and ip{x) ~ — ln(x/c) at x — +oo. The 
relation obtained ensures the renormalization invariance of as{q'^)- At low 
g"^, we obtain A^ = /i^ exp{— 47r/(6otts)}, which indicates the essentially non- 
perturbative character of three last terms in Eq. and these terms are 
absent in the usual perturbation theory. 

The acceptance of the cancellation mechanism for the nonphysical per- 
turbation theory singularity in Eq. (0) by the nonperturbative contributions 
leads to the necessity of supplementary definition of integral (H) near point 
k"^ = A"^. This problem can be reformulated as a problem of dividing pertur- 
bative and nonperturbative contributions in resulting in the introduction 
of some parameter feg ~ 1 GeV. This provides the absence of the pole at 
k'^ = A^ in both perturbative and nonperturbative parts Nonperturbative 
contributions in Eq. (0) decrease at infinity as the integral in Eq. (|^) 
converges and we can obtain 

K{A,ko) = Unn 

■^See also Ref. where the problem of perturbative and nonperturbative contributions 
to Qs is discussed and the definition of infrared finite regularized perturbative part of as 
is suggested. 



9 , 4:71 

g'An = — 
bo 
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Phenomenology gives the positive value of the gluon condensate K in the 
interval (0.32 Gel/)^ — (O.SSGeV)^ ||, |9|. If one takes the value of gluon 



condensate K^^'^ = 0.33GeV and consider formula ( [TT| ) as equation for A 
at different values of fco, the picture will be the following. At < = 
0.777 GeV^ there are no solutions in the interval A = 0.1 GeV — 0.434 Gel^. 
At /cq = ko there is a single solution A = "ilbMeV corresponding to ntg = 
rhg ~ 0.6 GeV. At fco > ko for A two solutions appear to whom two values mg 
correspond, one of them increases with the increase of ko and other decreases 
with increase of ko. This situation is illustrated in Fig. ^ where we used 
Eq. (|T0|) to connect the parameters A and 



m 
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It is seen from Eq. (^ that the pole singularities are situated at two points 
= and = —Tn?g- It corresponds to the two effective gluon masses, 
and TUg. Therefore, the physical meaning of the parameter mg is not the 
constituent gluon mass, but rather the mass of the exited state of the gluon. 

In a number of cases of the QCD calculations it is necessary to estimate 
integrals of the form 

F(g2) = j a,{e)f{e)dk, (12) 



where /(fc^) is some smooth function and integration includes the infrared 
region where the usual perturbation theory is inapplicable. In this cases 
the infrared matching scale fij can be introduced (A << /i/ << q) and 



the contribution to integral (0) from the region k > can be calculated 



perturbatively. To estimate the contribution of the region k < m Ref. 
the running coupling constant is assumed to be regular at zero and power 
terms are extracted using dimensional arguments. In this connection the 
nonperturbative parameters ap are introduced, 

p + 1 /■ 

«p(/^/) = — ^TT / dkasik)kP, (13) 
i 

which should be found from experiment. In the approach connected with 
infrared renormalons |31|, the analogous ambiguity arises when evaluating 



integral (p^. This ambiguity can be understood in the following way. Ex- 
panding cxsik"^) in the series of powers of as(g^) and integrating in k^ one 
obtains an asymptotic series with terms growing as factorial which can be 
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estimated by the finite number of terms up to power corrections. Tlie cliarac- 



teristic integral = ao{jj)/Ti (zero moment) lias been estimated [30, 32 
with the result aQ{2GeV) = 0.18 ± 0.02. This integral turned out to be 
not only stable with respect to the choice of the infrared regularization (fit- 
invariance) but also relative to different approximations of the high energy 
tail of the running coupling (number of active quarks, one- or two-loop ap- 
proximation). It is shown in Refs. |^ that the behaviour (|^) complies with 
integral estimations of Refs. Let us see whether the behaviour (|^) 



is compatible with these estimations. In this case we also encounter the 
problem of definition in expressions of form (|1^) and the necessity to fix the 



corresponding ambiguity. The following method seems to be highly conve- 
nient and universal, 

' li-:^77L(9.) • (14) 




(F)™ .'^0 de' (P)™ 

The finite arbitrariness described by the constant /i^ arises only if we integrate 
the logarithmic singularities (which are usually connected with physics). In 
the case of local integrability {n > 2m, here n is the dimension of integration 
space) the rule ( p^lf ) runs idle and at n < 2m it corresponds to analytic 



continuation in m. In accordance with the rule ([T^), the contribution of the 
singular term of Eq. (^ to the moments ([l3D is the following: 



Thus, at p < 1 the singular term of Eq. (|^ should be considered as a distri- 
bution and fit-invariance of the lowest moment aol^j) for different variants 
of infrared behaviour can probably indicate the absence of "freezing" the 
running coupling constant and universality of ambiguity fixing. In favour 
of "freezing" the coupling could indicate the equality of moments ap(/x/) for 
different p. For contribution of the last two terms of Eq. (|^ to zero moment 
one obtains: 



Aaoifii)/n = --—- --(c - if'arctgi^V^)- (16) 

Taking a ~ 0.4:2GeV, = 2GeV one finds rather small contribution of 
the singular term, A'^''3ao{^j)/TT = -(3a2)/(27r/i?) ~ -0.021. At A Ai 
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the second term of Eq. (|T6|) vanishes as — (16/9)(Ai — A)^/Af and it can be 
made small by a corresponding choice of A. For example, at A = 0.375GeV 
(in this case c = 1.3476) this term equals —0.0216. Having in mind the 
results of Refs. where the contribution of the first two terms of Eq. (|^) 
has been evaluated, we conclude that the running coupling constant Eq. (0) 
can be consistent with integral estimations in the infrared region. Note, that 
if we compare Eq. (5) and Eq. (12) of the second of Ref. [§, we find that 
" analytization" procedure in two- loop case leads, at large momentum, to 
the half nonperturbative contribution of one-loop case. It can point to the 
tendency of minimization of nonperturbative contributions in the ultraviolet 
region. 
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Figure 1: Momentum behaviour of the running couphng constant. Contin- 
uous curve corresponds to our Eq. (7); the curve (a) is one-loop Eq. (1); the 
curve (b) shows " analytized" Eq. (2); the curve (c) is "freezed" behaviour 
Eq. (3) where m is taken 1 GeV; the curve (d) corresponds to Richardson 
parameterization for the running couphng constant (Eq. (3) with = A). 
Here A = 300 MeV, bo = 9. 
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Figure 2: The dynamical gluon mass nig corresponding to the nonpertur- 
bative scale for fixed value of the gluon condensate X^/^ = 0.33 GeV. 
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